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Abstract 

We examine the CP-conserving (CPC) and CP-violating (CPV) effects of a 
general HZZ coupling through a study of the process H — > ZZ^ — > 
at the LHC. We construct asymmetries that directly probe these couplings. 
Further, we present complete analytical formulae for the angular distributions 
of the decay leptons and for some of the asymmetries. Using these we have 
been able to identify new observables which can provide enhanced sensitivity 
to the CPV HZZ coupling. We also explore probing CP violation through 
shapes of distributions in different kinematic variables, which can be used for 
Higgs bosons with tuh < 2 mz- 



1 Introduction 



The Standard Model (SM) has had unprecedented success in passing precision tests at the 
SLC, LEP, HERA and the Tevatron. However, the verification of the Higgs mechanism, 
which allows the generation of particle masses for fermions and electroweak (EW) gauge 
bosons without violating the gauge principle, is still lacking. The search for the Higgs 
boson and the study of its properties will be among the major tasks of the Large Hadron 
Collider (LHC), which will soon start operation, and of the International Linear Collider 
(ILC), which is under planning and consideration [1]. 

However, the instability of the Higgs boson mass to radiative corrections and the re- 
sulting fine tuning problem point towards the existence of physics beyond the SM (BSM) 
at the TeV scale. This BSM physics usually implies more Higgs bosons and may have 
implications for the properties of the Higgs boson(s). Hence, the determination of the 
Higgs boson quantum numbers and properties will be crucial to establish it as the SM 
Higgs boson [2] or to probe any new BSM physics. 

Furthermore, there is no real theoretical understanding of the relative magnitudes 
and phases of the different fermion mass parameters in the SM, even though we have 
an extremely successful description of all observed CP-violation (CPV) in terms of the 
Cabbibo-Kobayashi-Masakawa (CKM) matrix. Indeed, the CPV of the SM, observed 
only in the K -K and B -B systems to date, appears insufficient to explain the Baryon 
Asymmetry of the Universe (BAU) [3], and an additional source of CPV beyond that of the 
SM may be needed for a quantitative explanation. An extended Higgs sector together with 
CPV supersymmetry (SUSY) is one possible BSM option that may explain this BAU [4]. 
Thus it is clear that the knowledge of the properties of the Higgs sector and any possible 
CPV therein is of utmost importance in particle physics phenomenology at present [5,6]. 

The LHC will search for the SM Higgs boson in the entire mass range expected theoret- 
ically and still allowed experimentally [7,8], whereas precision profiling of the Higgs boson 
is expected to be one of the focal points at the ILC [9]. After discovery, the determination 
of the Higgs boson couplings, in particular those with a pair of electroweak gauge bosons 
(V = W/Z) and those with a pair of heavy fermions (/ = t/r), will be essential. In this 
study we focus on the HZZ coupling. 

The ILC, in both the e + e" and the 77 [10] options, and the LHC offer a wealth of 
possibilities for the exploration of the CP quantum numbers of the Higgs boson H [11]. 
At an e + e~ collider, the Z boson produced in the process e + e~ — > ZH is at high energies 
longitudinally polarised when produced in association with a CP-even Higgs boson and 
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transversely polarised in case of a CP-odd Higgs boson. The angular distribution of the 
Z boson therefore carries a footprint of the Higgs boson's CP properties [12-14]. Further- 
more, measurements of the threshold excitation curve can yield useful information on the 
spin and the parity of the Higgs boson and establish it to have spin and be even under 
parity transformation, hence J p = + , in a model-independent way [15,16]. Additionally, 
kinematic distributions of the final state particles in the process e + e~ — > ffH, produced 
via vector boson fusion or Higgsstrahlung, where / is a light fermion, with or without ini- 
tial beam polarisation, can be exploited to study the HZZ coupling, including CPV [13], 
[17-24]. Ref. [22] uses the optimal observable technique whereas Refs. [19,23,24] exploit 
the kinematical distributions to construct asymmetries that are directly proportional to 
different parts of a general CP-violating coupling. Associated production with top quarks 
e + e~ — > ttH may be used to extract CP information too [25,26]. 

Higgs decays may also be used effectively. The angular distributions of the Higgs 
decay products, either a pair of vector bosons or heavy fermions that further decay, can 
be exploited to gain information on the Higgs CP properties if it is a CP-eigenstate and 
the CP-mixing if it is CP violating [19], [27-30]. A detailed study of the Higgs spin and 
parity using the angular distributions of the final-state fermions in H — > ZZ — > leptons, 
above and below the ZZ threshold, was performed in [30]. The H — > // pair (/ = t/r) 
has the advantage of being equally sensitive to the CP-even and CP-odd part of the Higgs 
boson [31]. For Higgs bosons produced in association with heavy fermions, or Higgs decays 
to heavy fermions at an e + e~ collider, angular correlations and/or the polarisations of the 
heavy fermions may also be used [26,32,33]. 

An ILC operating in the 77 mode offers an attractive option not only for the CP- 
determination of the Higgs boson, but also for the measurement of a small CP-mixing in a 
state that is dominantly CP-even. Using linear and circular polarisation of the photons one 
can get a clear measure of the CP mixing [34]; further using a circular beam polarization, 
the almost mass degenerate CP-odd and CP-even Higgs bosons of the MSSM may be 
separated [35-39]. Interference effects in the process 77 — > H — > // (f — t/r) [40-44] can 
be used to determine the ffH and 77^ couplings for an H with indefinite CP parity. 

Hence, the e + e~ collider and its possible operation as a 77 collider offer some unique 
possibilities in the exploration of the CP quantum numbers of the Higgs boson. How- 
ever, the LHC is the next collider to come into operation. So we want to seek answers 
to these questions already at the LHC [45]. Here, the tt final state produced in the de- 
cay of an inclusively produced Higgs boson can provide knowledge of the CP nature of 
the ttH coupling through spin-spin correlations [46,47] whereas ttH production allows a 
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determination of the CP-even and CP-odd part of the // couplings with the Higgs boson 
separately [48,49]. The use of r polarisation in resonant t + t~ production at the LHC 
has also been recently investigated [50]. The HZZ coupling can be explored at the LHC 
in the Higgs decay into a Z boson pair which then decay each into a lepton pair, i.e. 
H -> ZZ& -> (e + £-)(£' + £'-) [30], [51-53]; above threshold, angular distributions have 
to be used while below threshold, the dependence on the virtual Z* boson's invariant 
mass may be exploited. Furthermore, this coupling (and the HWW coupling) can be 
studied in vector boson fusion [54-56], and a similar idea may be employed in H + 2 jet 
production [57,58] in gluon fusion (however, also see Ref.[59]). 

Most of the suggested measurements should be able to verify a scalar Higgs boson when 
the full luminosity of 300fb _1 is collected at the LHC (or even before), provided the Higgs 
boson is a CP eigenstate. For example, using the threshold behaviour it may be possible to 
rule out a pure pseudoscalar state with 100fb _1 in the SM [30]. However, a measurement 
of the CP mixing is much more difficult, and a combination of several different observables 
will be essential. 

In this paper we investigate CP mixing in the Higgs sector using the process, 
H _> ZZ& -> (£+£-) (£'+£'-). We extend the analysis of Ref. [30] to a Higgs boson of 
indefinite CP. Further, we extend the analysis of Ref. [53], where asymmetries were con- 
structed using angular distributions of the decay leptons, which directly probe the CP 
mixing. 

The paper is organised as follows. In section 2 we present the complete analyti- 
cal formulae for the angular distribution of the decay leptons produced in the process 
H -> ZZ& -> (e + £-)(£' + £'-), parameterising the HZZ vertex in a model-independent 
way, for a Higgs boson of indefinite CP. In section 3 we examine how this modified coupling 
changes the total number of H — > ZZ — > 4 lepton events seen at the LHC. In section 4 we 
then construct different observables that can be used to probe the CP nature of the Higgs 
boson and present the numerical results. In section 5, we propose an investigation of CP 
mixing using kinematical distributions of the decay leptons, and in section 6 we present 
our conclusions. 
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2 Model independent analysis of H — ► ZZ^ 

For our study of possible CPV in the Higgs sector we will examine the decay of a Higgs 
boson into two Z bosons with subsequent decay into two lepton pairs, 

H ^ZZ^^{hh){f 2 f 2 ). (1) 

To perform a model-independent analysis we examine the most general vertex including 
possible CPV for a spin-0 boson 1 coupling to two Z bosons with four-momenta q\ and q 2 , 
respectively. This can be written as 



HZZ ~ cos 9 



w 



(2) 



where p = q\ + q 2 and k = q\ — q 2 , 9w denotes the weak-mixing angle and e^ap is the 
totally antisymmetric tensor with €0123 = 1- As can be inferred from Eq. (2) the CP 
conserving tree-level Standard Model coupling is recovered for a = 1 and b = c = 0. 

The terms containing a and b are associated with the coupling of a CP-even Higgs 
boson to a pair of Z bosons, while that containing c is associated with that of a CP- 
odd Higgs boson. In general these parameters can be momentum-dependent form factors 
that may be generated from loops containing new heavy particles or equivalently from 
the integration over heavy degrees of freedom giving rise to higher dimensional operators. 
The form factors b and c may, in general, be complex. Since an overall phase will not 
affect the observables studied here, we are free to adopt the convention that a is real. This 
convention requires the assumption that the signal and background do not interfere, and 
indeed in our approximation where the Higgs boson is taken on-shell, this interference is 
exactly zero. Interference would be only manifest if the Higgs boson were taken off-shell 
and since the dominant signal contribution arises from on-shell Higgs bosons, we expect 
this interference to be small and neglect it. 

In principle, the vertex is valid at all orders in perturbation theory. Contributions to 
the HZZ vertex from loop corrections will not add any new tensor structures and will only 
alter the values of a, b and c. More generally, a, b and c are momentum dependent form 
factors obtained from integrating out the new physics at some large scale A. Since the 
momentum dependence will involve ratios of typical momenta in the process to the large 
scale A, we make the reasonable assumption that the scale dependence can be neglected 
and keep only the constant part. 



1 In fact, in order to be as general as possible one should allow for a general CP violating coupling with 
a "Higgs" particle of arbitrary spin, as in [30]. We keep this for future work. 
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Non-vanishing values for either Qm(b) or 9m(c) destroy the hermiticity of the effective 
theory. Such couplings can be envisaged when going beyond the Born approximation, 
where they arise from final state interactions, or, in other words out of absorptive parts 
of the higher order diagrams, presumably mediated by new physics. Further, a, dte(b) 
and ^sm(c) are even under T, while 5sm(b) and 9fe(c) are odd, where T stands for the 
pseudo-time reversal transformation, which reverses particle momenta and spins but does 
not interchange initial and final states. It is the CPT odd coefficients that are related to 
the presence of absorptive parts in the amplitude [60]. In most CPV extensions of the 
SM one has \a\ 3> \b\, |c|, so most of the observables used to study the HZZ vertex are 
dominated by the first term in the vertex Eq. (2); in order to probe the last, the CP-odd 
term, it is most advantageous to construct asymmetries which vanish as CP is restored. 

CP violation will be realized if at least one of the CP-even terms is present (i.e. either 
a^O and/or b ^ 0) and c is non-zero. In the following we keep the three coefficients non- 
zero in our analytical work, where appropriate. However, in the numerical presentation 
of most of our results we will take b = for simplicity, keeping non-zero b only where 
essential. Further, we make the justified approximation to neglect the possible momentum 
dependence of the form factors. 

Notice that neither qx^^zz nor QzvVhzz are zero > i- e - the Ward identities are violated. 
This is due to the breaking of electroweak symmetry and is already the case for the SM 
vertex. Some studies, e.g. Refs. [24,53], explicitly construct the extra terms such that they 
satisfy such Ward identities individually, for example, by taking a CP-even term of the 
form q 1 ■ q^g^v — qi^qiv Strictly speaking, this is not necessary as long as any additional 
terms vanish in the limit vrtz — > 0. Furthermore, since one must separately include the SM 
g^ v coupling and the new CP-even contribution (with independent coefficients), one may 
always reproduce our choice of the vertex with a suitable redefinition of the coefficients. 

Our vertex differs from the vertex of Refs. [15,30] only in the choice of the normalisation 
of the coefficients (to make them dimensionless) . The normalisation of the coefficients (and 
the overall normalisation) also differs from Refs. [51,56], where uih was used in contrast 
to our m Z - Additionally, Refs. [51,56] use the momenta of the Z-bosons to define the last 
term (i.e. ~ q^q^ ) in contrast to our ~ p a k fi . However, this last difference is for this process 
only a factor of —2 since the additional terms are removed by the asymmetric property 
of the tensor. Finally, Ref. [24] differs in the choice of the last term (again ~ qfq^ ) and 
rearranges the contributions of the first two terms, as discussed in the preceding paragraph. 
For 6 = and light lepton final states, all these vertices are the same, modulo momentum 
independent normalisations of the coefficients. 
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From the above discussion it is clear that the total decay rate of Eq. (1), which is 
CP-even and T even, can only probe a, dte(b) and the absolute values of b and c. In 
order to probe the other non-standard parts of the HZZ coupling, in particular in order 
to probe CP-violation, one must construct observables that are odd under CP and/or 
T. These observables give rise to various azimuthal and polar asymmetries and will make 
their presence felt through rates which are integrated over a partial (non-symmetric) phase 
space. Thus one may probe 9fte(6), £sm(b), 3?e(c) and Qm(c) either by using the shapes 
of various kinematical distributions or by constructing observables which are obtained 
using partially integrated cross sections [19,23,24] 2 . We will use the latter to construct 
asymmetries which receive contributions from non-standard couplings and which vanish in 
the tree-level SM. These are related to simple counting experiments, recording the number 
of events in well defined regions of the phase space. It may also be noted that results 
obtained using these asymmetries are less sensitive to the effect of radiative corrections to 
the production [61] and decay [62,63] of the Higgs boson. 

In order to find observables which project out the various non-standard couplings in 
Eq. (2) it is instructive to have an analytical formula for the differential distribution of the 
Higgs decay to off-shell Z bosons with subsequent decay into fermion pairs with respect to 
the various scattering angles. We denote the polar angles of the fermions fi, f 2 in the rest 
frame of the parent Z bosons by Q\ and 6> 2 , and the azimuthal angle between the planes 
formed from the fermion pairs in the Higgs rest frame by [see Fig. 1]. Also note that 
there can be no angular correlations (at tree-level) between the initial and final states (i.e. 
between the beam-direction and the final state leptons) as long as the Higgs has zero spin. 



Introducing the notation cg i = cos 6>j, s$ i = sin 6>j {i = 1,2), = cos0, etc. the 
tree-level differential decay rate for distinguishable fermions can be cast into the form 



dce 1 dcg 2 d(j) 



~ a 



,2 „2 



■ S V#, - TT-^ei^e^ + [(1 + cgj(l + c 2 02 ) + s 2 ei s 2 62 c 2< p] 



2r?i^2 ( 1 
I se 1 se 2 c <t> c 6»i c 6> 2 

la \ la 



+ 1 1 2 X S 9i S d 2 
! a 



+ C 



2% 



Ax 1 



2 In fact, Ref. [24] constructed systematically the whole set of asymmetries which probe different parts 
of the anomalous couplings. 



6 



Figure 1: The definition of the polar angles 6i (i = 1, 2) and the azimuthal angle <fi for the 
sequential decay H — > ZZ^*> — > (/i/i) ^2/2) ■ 
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(3) 



where 2 = mim 2 /m 2 z with mi, m 2 the virtualities of the Z bosons (q 2 = m 2 ). Furthermore, 
we have introduced the notation 7 a = 7172(1 + P1P2) and % = 7i72(A + P2) m terms of 



(4) 



the Lorentz boost factors of the Z bosons, 7* = l/y 1 — /3f , and the velocities 

A = f = l,2, 

where ^ are the Z boson energies in the Higgs rest frame and 



(mi + m 2 ) 2 *_(m l — m 2 ) 



m 



H 



ni 



H 



1/2 



The rji are given in terms of the weak vector and axial couplings Vf v a/\, 

■ 2rf " L with v fi =Tl-2Q fi sm 2 d w , a fi = T% 



2 1 2 ' 
vi + ai 

Ji Ji 



(5) 



(6) 



Here T 3 . denotes the third component of the weak isospin and Qf i the electric charge of 
the fermion f iy in our case e" or 



3 Sensitivity of the total production to new couplings 

As discussed in section 2, one may use the total decay rate of the process in Eq. (1) to test 
possible deviations from the SM in the Higgs to ZZ coupling. At the LHC the dominant 
Higgs production process is given by gluon-gluon fusion, 

gg^H^ZZV^ifJMhh), (7) 

with / = e or /i. The width for the process H — > ZZ^ — > (/i/i) (fzfi) is given by, 
T(H -> ZZ« -> (A/0 (/ 2 / 2 )) = 

vr 2 io 1 io 2 [K - ™!) 2 + m|r|] [(ml - ml)- + m|r|] MZ ' W 

where the width for the Higgs decay to two Z bosons 3 of virtualities mi and m 2 is, 



h^zz — 



16V27T 



(3 < a 



,2 



„ 9 I2m 2 ml 
/r + 



m 



H 



+ \b\ 2 ^f ^ + \c\ 2 x 2 8(3 2 
m% 4 



2 -\ 

+aKe(6) /3 2 ^ (3 2 + 4m 2 m 2 /m|, J (9) 



and T z ^jj t is the width for the decay of a Z boson to a fermion pair, /j/j, as given in the 
SM, 

r ^ = f^^(^ + 4)- ( 10 ) 

As expected, the CPT-even total rate cannot directly test CPV (since there is no 
interference between the CP-even and CP-odd terms), but it is sensitive to possible non- 
SM coupling effects in 3t.e(b) and the absolute values of b and c. Furthermore, Eq. (9) 
shows that the linear rise in (5 just below the threshold is typical [15] of the SM Higgs 
boson 4 . 

The Tevatron is in principle also sensitive to the process of Eq. (7) for sufficiently 
high Higgs boson masses. Indeed, preliminary Tevatron results [65] indicate that a sig- 
nal for a Higgs boson of 150 GeV would have been seen (with 95% confidence) if the 



3 For the on-shcll decay H ZZ, sec Rcf. [64]. 

4 This observation is valid for all spins, with one minor caveat: the spin-2 case can also have a term 
which presents a linear rise in but this can be excluded by angular correlations, see Ref. [15]. 
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observed(expected) DO-CDF combined total cross-section were enhanced by a factor of 
2.4(3.3). However, this result is dominated by the decay H — > W + W~; the H — > ZZ de- 
cay is suppressed relative to W + W~ by around a factor of 10 for a 150 GeV Higgs boson, so 
an enhancement of the HZZ vertex from additional couplings would need to be very large 
indeed to be seen by the Tevatron. Since we are here investigating the HZZ coupling, we 
make the assumption that the other decay channels are unaffected and that any change 
originates from the HZZ coupling alone. For lower Higgs masses the HZZ coupling can 
also play a role in the production of the Higgs via the channel qq — > Z* — > ZH. However, 
as can be seen from Ref. [66], with current data, the Tevatron would be sensitive to this 
production mode only if the cross-section were enhanced by a factor of ~ 30 — 90 compared 
to the SM and thus the nonobservation of this channel in the current data only puts very 
weak constraints on the magnitude of these couplings. 

To estimate the sensitivity of the LHC to deviations from the SM coupling, we refer 
to the ATLAS study for the process of Eq. (7) at m H = 150 GeV and 200 GeV [7,45]. 
In this study, four leptons were selected using the standard electron and muon identi- 
fication criteria. Events were required to have two leptons with pr > 20 GeV and two 
additional leptons with p T > 7 GeV, with rapidity \rj\ < 2.5 for all four. The signal and 
background were compared in a small mass window around the Higgs boson mass, and a 
lepton identification and reconstruction efficiency was applied. 

For the tjih = 150 GeV analysis, one lepton pair was required to have an invariant mass 
within 10 GeV of mz while the other pair was required to have an invariant mass above 
30 GeV. Additionally, isolation and impact parameter cuts were used to further remove irre- 
ducible backgrounds. For the 200 GeV analysis, the continuum ZZ background was further 
removed by requiring the pt of the hardest Z-boson to be greater than m#/3 ~ 66.6 GeV 
(see Refs. [7,45] for further details). 

Note that this ATLAS study was performed at tree-level with no K-factors. Higher 
order corrections to the production process could alter the cross section by up to a factor 
two [61]. The higher order electroweak corrections to the Higgs decays into W/Z bosons 
have been calculated in Ref. [62] in the narrow width approximation. Ref. [63] presents 
the complete 0(a) corrections to the general H — > 4/ processes, including off-shell gauge 
bosons which are important for our study. The corrections have been shown to change 
the partial width by up to 5% for the Higgs boson masses we consider in this paper. Our 
analysis, which uses the results of the ATLAS study, strictly speaking is only valid at 
tree-level, despite the all-orders validity of the HZZ coupling (see section 2). 

After these cuts, the study found for a 150 GeV Higgs boson and an integrated lumi- 
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nosity of 100 fb^ 1 , 67.6 signal events with a background of 8.92 events. The corresponding 
signal and background events for a 200 GeV Higgs boson were 54 and 7, for an integrated 
luminosity of 30 fb" 1 . Altered HZZ couplings will enhance (or decrease) the number of 
signal events, while leaving the number of background events fixed. However, the size of 
this enhancement (or reduction) is model-dependent. Although the change in the width for 
H — > ZZ* — > 4/ is clear from Eq. (9), the branching ratio depends on how the other Higgs 
decay channels are affected by the new physics. As mentioned above, we here make the 
assumption that only the HZZ vertex deviates from that of the SM. If this were not the 
case, and, for example, the HWW coupling was similarly enhanced, then any enhancement 
of the H — > ZZ branching ratio would be watered down. Furthermore, we assume the 
Higgs production proceeds as according to the SM, since the dominant production mode 
contains no HZZ coupling, but one should be aware that CPV effects in other vertices 
may alter the Higgs production rate (see e.g. Ref. [67]). Finally, we assume that the rate 
calculated with the general HZZ coupling Eq. (2) will be reduced by experimental cuts 
in the same way as the SM rate. Only electron and muon final states are considered, and 
we scale up the number of signal and background events to correspond to an integrated 
luminosity of 300 fb _1 . 

We then calculate the total number of signal events N s that we expect from the 
new coupling and compare the expected change (with respect to the SM) with the 
possible statistical fluctuations of the SM signal and backgrounds. The significance of 
this deviation from the SM expectation (in units of one standard deviation) is then 
(Ns — Af M )/ a/A| m + Nb, where N$ M is the number of signal events expected in the 
SM and Nb is the number of background events. This quantity, for m# = 150 GeV and 
200 GeV, is plotted in Fig. 2, where we have scanned over values of the couplings a and 
|c| (the total rate is independent of the phase of c). For simplicity, we have set b = (for 
6 = 0, one can see that Eq. (9) is symmetric in a allowing us to restrict the plot to positive 
values). As can be inferred from Fig. 2 in the white region we can not distinguish the 
corresponding a, c values from the SM 1, c = at a significance more than 3 a. 

Large values of |c|, however, together with the SM value of a = 1 are easily identified at 
the LHC. For example, the scenario a = c = 1 is excluded with around 5 a significance for 
m H = 150 GeV and over 20 a significance at m H = 200 GeV. However, since |c| arises from 
new physics one would expect its value to be suppressed by the size of the new physics scale, 
and therefore be rather small. For a = 1 (the SM value) we find that this measurement 
provides 3 a evidence of non-zero c only if c > 0.75 or c > 0.32, for m H = 150 GeV and 
200 GeV, respectively. Furthermore, since both a 2 and |c| 2 contribute to the total rate, we 
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Figure 2: The number of standard deviations from the SM which can be obtained in the 
process gg — > H — > — > 4 leptons, as a scan over the (a, |c|) plane. The Higgs mass 
has been chosen to be 150 GeV (left) and 200 GeV (right). The white region is where the 
deviation from the SM is less than 3 a; in the light blue/light grey region the deviation 
is between 3 a and 5 a; while for the dark blue/dark grey region the deviation is greater 
than 5 a. 

cannot distinguish whether or not any deviation is originating from non-standard values of 
a or |c|, and even if the SM total rate is confirmed, one cannot definitively say that a and c 
take their SM values since an enhancement in |c| may be compensated by a reduction in a. 
Also, a non-zero value of b could provoke a similar effect. Indeed, the total rate is not even 
reliable in distinguishing a CP-even eigenstate from a CP-odd one. Instead, to provide 
a definitive measurement of CP violation in this coupling, one must explore asymmetries 
which probe the interference of the CP-even and CP-odd contributions directly. 

4 Asymmetries as a probe of CP-violation 

As stated above, apart from the terms proportional to a and dte(b), all other contributions 
to the vertex Eq. (2) are odd under CP and/or T transformations, and their presence 
implies violations of the corresponding symmetries in the interaction. We exploit this 
by constructing observables from the 3-momenta of the initial and final state particles 
with the same transformation property under the discrete symmetries as one of these 
non-SM couplings. The expectation value of the sign of such a variable will directly 
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probe the corresponding coupling coefficient [24]. 5 The asymmetry will be proportional 
to the probed coupling and therefore non-zero only if the corresponding non-SM coupling 
is present. Furthermore, since these asymmetries are exactly zero for all backgrounds 
(we neglect interference effects), backgrounds cannot contribute to the asymmetry, except 
through fluctuations, and it is therefore possible to use less stringent cuts on the signal. 

In this section we present various observables and their asymmetries which allow one to 
probe the real and imaginary parts of the form factors b and c, the latter being indicative 
of CP violation for simultaneously non-zero a and/or b values. 

1. An observable to probe Ssm(c): We consider the observable 



subscripts Z and H denote that the corresponding 3-vector is taken in the Z boson or 
Higgs boson rest frame, respectively. This observable is CP odd and T even and thus 
probes the non-SM coupling with the same transformation properties, i.e. 3m(c). With 
the above angular definitions we have 



We can calculate the resulting asymmetry by integrating Eq. (3) over the angles with an 
appropriate weighting. Although Eq. (3) is only valid for distinguishable fermions, we may 
include fermions of the same flavour, e.g. (e~e + )(e~e + ), and distinguish the fermions by the 
requirement that the first pair reconstruct the Z-boson mass. In general, the contribution 
from the same final state with the antiparticles switched would contain two off-shell Z- 
bosons and may be neglected. However, one should also note that this observable requires 
one to distinguish between fermions and ant i- fermions. 

The angular distribution of Eq. (3) contains several terms linear in cos#i. However, 
most of these terms are removed by integration over the angles 9 2 and 0, leaving only 
one term proportional to aSm(c). So only a non-zero value of ^sm(c) gives rise to this 
forward-backward asymmetry and hence provides a definitive signal of CP violation in the 
HZZ vertex. This is demonstrated in Fig. 3, which shows the dependence on cos#i for 
pure CP-even, pure CP-odd and CP- violating interactions 6 . 

5 This statement is true strictly when only the linear terms in the anomalous HZZ coupling are kept. 
Potentially, the asymmetries may also contain combinations of more than one (small) anomalous couplings 
which will have the same discrete symmetry transformation properties. In that case the asymmetry will 
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Here p iy i = 1, . . .4 are the 3-momenta of the leptons (in the order /1/1/2/2), and the 



0\ = cos 61 . 



(12) 
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Figure 3: The normalized differential width for H — > — > (/i/i) (72/2) an d 
m-H = 200 GeV with respect to the cosine of the fermion /i's polar angle #1. The solid 
(black) curve shows the SM case (a = 1, b = c = 0) while the dashed (blue) curve is for 
a pure CP-odd state (a = b = 0, c = i). The dot-dashed (red) curve is for a state with 
a CP violating coupling (a = 1, b — 0, c = i). One can clearly see an asymmetry about 
cos 6*i = for the CP violating case. 

To quantify the effect we define an asymmetry by, 



This asymmetry, which is the expectation value of the sign of cos^i (Eq. 12) and which is 
CP-odd and T even, directly probes Qm(c) which is also CP-odd and T even. Integrating 
Eq. (3), the asymmetry A± can be written as 



where f is related to the decay width H — > ZZ^> — > (/1/1) (^2/2), c.f. Eqs. (8,9), and is 
given by 



and the integral is over the virtualities, weighted with the Breit-Wigner form of the Z-boson 

be a direct probe of that particular combination of the non-SM couplings. 

6 This figure differs from the corresponding figure in Ref. [53] for the CP violating coupling due to the 
different conventions. The corresponding curve for the mixed CP state in Ref. [53] is reproduced with our 

current conventions if a = 1, b = 0, c = —i/2. 



r(cos#i > 0) - r(cos#i < 0) 
r(cos#i > 0) + r(cos#i < 0) ' 



(13) 




(14) 




(15) 
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propagators, 
d 2 V ■ ■ ■ 



r m H r 
I dm\ 
Jo Jo 



[ntH-m,!] 



dm 



rrin 



[(m\ - mlf + m|r|] [(ml - m\f + m\Y\\ 



(16) 

This asymmetry is calculated at tree-level. Higher order electroweak corrections to the 
decay H — > ZZ — > 4 leptons are of the order 5-10% for angular distributions [62,63]. 
One might worry that these corrections could feed into the asymmetry and swamp the 
signal. However, unless the corrections introduce some new effect (and are thus in some 
sense "leading order"), one expects their contribution to CP violation to be of a similar 
proportion as those at tree-level, so they would provide a correction to Eq. (14) of 5-10%, 
and not significantly alter our results. 

Fig. 4 shows the values of A\ for a Higgs mass of 150 and 200 GeV, respectively, as 
a function of the ratio 3m(c)/a and where we have set b = for simplicity. The value 
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Figure 4: The asymmetry A\ given by Eq. (14) as a function of the ratio 3m(c)/a, for 
a Higgs boson of mass 150 GeV (left) and 200 GeV (right). We chose 6 = 0. The inserts 
show the same quantities for a larger range of ^sm(c)/a. 



Qm(c)/a = corresponds to the purely scalar state and Qm(c)/a — > oo to the purely 
CP-odd case. It is clear from Eq. (14) that A± is sensitive only to the relative size of 
the couplings since any overall factor will cancel in the ratio, c.f. Eq. (13). We find that 
the asymmetry is maximal for Qm(c)/a ~ 1.5(0.7) with a value of about 0.067(0.077) 
for mn — 150(200) GeV. The smallness of this asymmetry arises from the fact that it is 
proportional to the coupling r\\ = 2v\a\j(v\ + af) which is equal to approximately 0.149 
for e,fi final states, c.f. Eq. (14). 

In order to estimate whether this asymmetry can be measured at the LHC, we calculate 
the significance with which a particular CP violating coupling would manifest. To do this, 
we must take into account the backgrounds to the signal process, which will contaminate 
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the asymmetry in two ways. Firstly, despite being CP-conserving the backgrounds may 
contribute to the numerator of the asymmetry via statistical fluctuations (e.g. the back- 
ground events with 0\ > may fluctuate upwards while those with 0\ < may fluctuate 
downwards and vice versa). Secondly, they will directly contribute to the denominator of 
the asymmetry. 

Consequently, the measured asymmetry will be given by, 

iVT ym N s 
Al ~ N s + N B - Al N s + N B > [U) 

where iV| sym is the asymmetry in the number of events in the two hemispheres, and Ai is 
the perfect theoretical asymmetry given in Eq. (13). 

The statistical fluctuation in an asymmetry calculated using a total number of events 
iV = N B + N s , even when N B and N s are expected to be symmetric, is 1/y/N. Hence, the 
significance of the expected asymmetry, S, in units of this statistical fluctuation is given 
by 

S = A™ cas ^fN = -L = A l —L. (18) 

In order to calculate this, we need to know the number of signal and background events 
expected at the LHC. However, in this case, since the contamination of the significance 
from the background is rather minimal, we choose to use the event sample before the 
detailed cuts to remove backgrounds, but after the initial selection cuts. For 150 GeV we 
take the number of signal and background events before applying the additional isolation 
and impact parameter cuts to remove the irreducible backgrounds, and for 150 GeV we do 
not apply the final pr cut on the hardest Z-boson (see Refs.[7,45]). 

Then, according to Refs.[7,45], for a m H = 150 GeV SM Higgs boson, we have a signal 
cross-section of 5.53 fb, with an overall lepton efficiency of 0.7625. Assuming an integrated 
luminosity of 300 fb -1 this gives 1265 signal events. For m H = 200 GeV, the corresponding 
signal is 1340 events. The number of signal events for the CP violating case is then 
obtained by multiplying the number of SM events by the ratio of CP violating to SM 
branching ratios. In the CP-violating case we always assume the SM value for the CP-even 
coefficient, a = 1. For simplicity we assume the charge of the particles to be unambiguously 
determined, and pair the leptons by requiring at least one pair to reconstruct the Z boson 
mass. The number of background events before cuts has been derived correspondingly 
from the study Refs.[7,45] and amounts to 1031(740) events for m B = 150(200) GeV. 

The significances are shown in Figs. 5 for itlh = 150 and 200 GeV, respectively, as a 
function of Qtt^c) with a = 1 and 6 = 0. As can be inferred from the figures the maximum 
of the curves is slightly shifted to higher values of 9m(c)/o compared to the corresponding 
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Figure 5: The significances corresponding to the asymmetry A± as a function of Qro^c), 
for a Higgs boson of mass 150 GeV (left) and 200 GeV (right). We chose the CP-even 
coupling coefficient a = 1 and b = 0. The inserts show the same quantities for a larger 
range of 3m(c). 



Figs. 4. This is due to the increasing Higgs decay rate with rising pseudoscalar coupling. 
The curves show that, even in a best case scenario, the significance is always < 3.5 a. This 
asymmetry may provide only evidence for CP violation (i.e. a greater than 3 a deviation 
from the SM) if ^m(c) > 1.9(0.7) for m H = 150(200) GeV. 

However, since one does not need to distinguish f 2 and f 2 one could also consider using 
jets instead of muons, i.e. H — > ZZ — > l + l~jj, to increase the statistics. If we use the bb 
final state, one can benefit from the increase by a factor ~ 4.5 in the branching ratio of 
the Z boson into a bb pair relative to the branching ratio into a lepton pair. As a matter 
of fact a study by ATLAS [68] shows that for a Higgs boson mass of 150 GeV with 30 fb _1 
it is possible to have a Higgs signal with a significance of 2.7a in this channel. So indeed 
one can foresee the use of this channel to add to the sensitivity. 

2. Observables which probe 5?e(c) and/or 3?e(6*c): We have constructed several 
observables which allow one to probe Ke(c). For this we need an observable which is CP 
odd and T odd. One possible observable is given by 

\P2Z -PizWpah x p 3 h\ 
which in terms of the scattering angles reads 

2 = - sin sin 0i . (20) 

(Since sin#i is always positive, one could equivalently use sin0 as the observable and 
obtain the same results.) By comparing this angular dependence with the differential 
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angular decay width given in Eq. (3), one can see that the corresponding asymmetry 
should pick up the third term ~ 771772 of the contribution multiplied with c$le(c) and the 
second term of the contribution multiplied with 3?e(6*c) and which also contains 771772- And 
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Figure 6: The asymmetry A 2 given by Eq. (21) as a function of the ratio 9?e(c)/a, for a 
Higgs boson of mass 150 GeV (left) and 200 GeV (right). We chose 6 = 0. The inserts 
show the same quantities for a larger range of Ke(c)/a. 



indeed we find for this asymmetry 

T(0 2 > 0) - T(0 2 < 0) 



A 2 



r(o 2 >o) + r(o 2 <o) 

= 1 J d 2 V Vir]2 x lb [a^e(c) la + ^e(b*c) X1 l] . 



(21) 



By construction, for b = or to linear order in the anomalous couplings, it is proportional 
to 9fe(c) as expected. This asymmetry is plotted in Figs. 6 as a function of 9te(c)/a 
for m H = 150 and 200 GeV, respectively. Since the form factors b, c are expected to 
be small we do not expect terms of second order in these coefficients to have a large 
impact, so here and in the following we set b = 0. Indeed, for the asymmetry A 2 with 
$te(b*c) ~ 3?e(c) 2 < 0.5 the change in the asymmetry due to neglecting b is < 30%. 
Figs. 6 show that this asymmetry is very small, with values below about ~ 0.011, which 
is principally due to the proportionality to the small quantity 771772 in Eq. (21). The 
significances for the asymmetry A 2 are shown in Figs. 7 for the two Higgs boson mass 
values. With values below about 0.55 they are far too small to provide evidence for CP- 
violation due to non-zero 3?e(c). Furthermore, in this case one cannot exploit the decay of 
Higgs bosons to jets since one must also distinguish p$H and pm- 

The smallness of the asymmetries A\ and A 2 are directly due to their proportionality 
to the factors 771, 772 . Looking at Eq. (3), one sees that this is true for all terms proportional 
to a5m(c), so not much can be done to improve on A\. However, this is not the case 
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Figure 7: The significances corresponding to the asymmetry A 2 as a function of Ke(c), for 
a Higgs boson of mass 150 GeV (left) and 200 GeV (right). We chose the other coupling 
coefficients a = 1 and 6 = 0. The inserts show the same quantities for a larger range of 
3te(c). 



for terms proportional to a3?e(c). So we may take our cue from the explicit analytical 
expression to construct new observables for which the asymmetry will not have these 
suppression factors. One such observable is given in terms of the angles by 



3 = cos Q\ sin 6 2 cos 9 2 sin . 



(22) 



O3 can be rewritten using the definition of 0\, c.f. Eq. (11), in terms of the four three- 
vectors, 



where 



3a 
3b 



o 3 = o 1 o 3a o 3b , 



(PiZ ~ P3Z) ■ (PlH X p 2H ) 

Waz - PszWpih x P2H\ 
(Pzz - Paz) • (pih + P2h) 



(23) 



(24) 



\P3Z - PazWpih + P2H\ 

In order to exploit this observable, we have to discriminate between all four leptons. For 
the asymmetry A3, 

T(0 3 > 0) - T(0 3 < 0) 



A 3 



r(o 3 >o) + r(o 3 <o)' 



we find analytically 



A 3 = I J d 2 V (^) [a^e(c) la + ^e(b*c) X1 2 b ] 



(25) 



(26) 
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Note that it no longer contains the suppression factors rji, r] 2 and for b = it probes 
the real part of the form factor c. By comparing the angular structure of O3 with the 
differential angular distribution Eq. (3), one sees that the asymmetry A3 picks up the first 
term in the contribution proportional to aKe(c) and the first term in the one proportional 
to 3?e(6*c). A non-zero value of A3 is hence an unambiguous sign of CP- violation. 




Re(c)/a Re(c)/a 

Figure 8: The asymmetry A3 given by Eq. (26) as a function of the ratio 5fte(c)/a, for a 
Higgs boson of mass 150 GeV (left) and 200 GeV (right). We chose 6 = 0. The inserts 
show the same quantities for a larger range of §fte(c)/a. 



Figs. 8 show the asymmetry A3 for m H = 150 and 200 GeV, respectively, where we 
have taken b = for simplicity. With values of < 0.09 they are about a factor 10 larger 
than those of Ai- The corresponding significances which should be achievable at the LHC 
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Figure 9: The significances corresponding to the asymmetry A3 as a function of 3?e(c), for 
a Higgs boson of mass 150 GeV (left) and 200 GeV (right). We chose the other coupling 
coefficients a = 1 and 6 = 0. The inserts show the same quantities for a larger range of 
Ue(c). 



for this asymmetry are shown for m# = 150 and 200 GeV in Figs. 9. They are maximal at 
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Ue(c) « 3(1) for m H = 150(200) GeV. For a 150 (200) GeV Higgs boson this asymmetry 
would provide evidence for CP-violation for dte(c) > 1.25(0.6) though discovery (a 5er 
deviation) is still out of reach. 

One should note, however, that a zero value for this asymmetry does not imply the 
absence of CP-violation, since for b ^ it could also happen that the contributions pro- 
portional to aKe(c) and Ke(6*c) cancel and mimic CP-conservation. In order to unambigu- 
ously show CP-violation in the HZZ coupling we hence need an additional observable to 
determine the two unknowns 3?e(c) and Ke(6*c). Such additional observables are presented 
in the following. 

An observable, which probes 9fte(c) alone, is given by 

Q ^ = [{p3H x pah) ■ Pih][(P3H * Pah) • (pin x P2h)} ^ 

' \P3H +PAH\ 2 \PlH +P2H\\P3Z ~P4z\ 2 \PlZ -P2z\ 2 /^-Q ' 

In terms of the angles it reads 

4 = sin 2 6\ sin 2 62 sin <\> cos . (28) 

(Again, since sin 2 9i j2 are always positive, this is equivalent to using an observable sin 20.) 
This coupling structure appears in the decay width only in the contribution which is 
proportional to aKe(c), c.f. Eq. (3), so that we can expect the corresponding asymmetry to 
probe CP-violation due to simultaneous non-vanishing form factors a and c unambiguously. 
Indeed the asymmetry is given by 

r(o 4 > 0) - r(o 4 < 0) 
4 ~ r(o 4 > 0) + r(o 4 < 0) 

= 1 J d 2 V ^pjaMe(c)x'y b . (29) 

Furthermore, as can be inferred from Figs. 10, which show the asymmetry for m H = 150 
and 200 GeV as a function of 5fte(c)/a, the asymmetries are larger than those of A2 and 
A3, with maximal values of up to ~ 0.11. The significances which may be achieved at the 
LHC are shown in Figs. 11. They reach values of up to almost 5 for m H = 150, 200 GeV 
so that this observable may probe CP-violation in an unambiguous way at the LHC for 
sufficiently large values of 9te(c). As can be inferred from Figs. 11, for a 150 GeV Higgs 
boson evidence for a non-zero 9te(c) is possible for 3?e(c) > 1, while for a 200 GeV Higgs 
boson evidence is already possible for 9fe(c) > 0.4. 

Alternatively one may use a combination of O3 and 4 to test CP-violation due to non- 
vanishing aKe(c) and/or Ke(6*c). For example, in terms of the angles a possible observable 
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Figure 10: The asymmetry At given by Eq. (29) as a function of the ratio Ke(c)/a, for 
a Higgs boson of mass 150 GeV (left) and 200 GeV (right). The inserts show the same 
quantities for a larger range of Ke(c)/a. 
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Figure 11: The significances corresponding to the asymmetry A4 as a function of 5te(c), 
for a Higgs boson of mass 150 GeV (left) and 200 GeV (right). We chose the other coupling 
coefficients a = 1 and 6 = 0. The inserts show the same quantities for a larger range of 
$fte(c). 



O5 is given by 



O5 = sin 6\ sin 62 sin 0[sin 6\ sin 62 cos <j) — cos 6\ cos 6*2] 



and can be constructed from the three-vectors by 

[(P4H X p 3H ) ■ Pih] [(PlZ - P2Z) ■ P3Z] 



The related asymmetry 



\P3H + P4H\\P3Z -P4z\ 2 \PlZ ~P2z\ 2 /8 



A 5 



T(Q 5 > 0) - r(o 5 < 0) 

r(o 5 > 0) + r(o 5 < 0) 



(30) 



(31) 



(32) 
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is shown in Figs. 12 for m# = 150, 200 GeV and yields the largest values among the asym- 
metries discussed so far, up to ~ 0.15. In Figs. 13 we show the corresponding significances. 
We see that for a 150 (200) GeV Higgs boson, this asymmetry would provide evidence for 
CP-violation for &e(c) > 0.66 (0.25) and discovery of CP-violation for 3te(c) > 1.28(0.52). 




Figure 12: The asymmetry A$ given by Eq. (32) as a function of the ratio 5fte(c)/a, for 
a Higgs boson of mass 150 GeV (left) and 200 GeV (right). We chose 6 = 0. The inserts 
show the same quantities for a larger range of 9?e(c)/a. 
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Figure 13: The significances corresponding to the asymmetry A5 as a function of dte(c), 
for a Higgs boson of mass 150 GeV (left) and 200 GeV (right). We chose the other coupling 
coefficients a = 1 and 6 = 0. The inserts show the same quantities for a larger range of 
£e(c). 



3. An observable which probes Sm(6): For completeness, we also present an ob- 
servable that probes the imaginary part of the CP-even form factor b. It is given by the 
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following combination of three-vectors 

[{piz - P2z) ■ (Pm + Pah)} [{p3H x Pah) • Pih] 



\piz — P2z\ 2 \psh + Pah \ 2 \pzz ~Paz\/^ 
= sin 6\ cos 6\ sin 6 2 sin <p . 



(33) 



And the asymmetry reads analytically 

r(o 6 > o) - r(o 6 < o) 



A = 



r(o 6 >o) + r(o 6 <o) 



i 
r 



d 2 V -7]2 a^sm(b) x'j 2 
8 



(34) 



Figs. 14 and 15 show the corresponding asymmetries and significances. Notice that once 
again, the asymmetry is proportional to the small factor i] 2 and is therefore rather small, 
i.e. < 0.025. Correspondingly this observable does not provide a good significance (only 
reaching values of about 1), so that the extraction of ^m(b) from this observable does not 
seem to be feasible at the LHC. Unfortunately, since this small factor is present in all the 
relevant terms in Eq. (3), all asymmetries that one can construct to probe this coefficient 
will be similarly small. 
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Figure 14: The asymmetry Aq given by Eq. (34) as a function of the ratio ^sm(b)/a, for 
a Higgs boson of mass 150 GeV (left) and 200 GeV (right). We chose c = 0. The inserts 
show the same quantities for a larger range of ^sm(b)/a. 



Refs. [28,29] also consider reweighting observables with the product of the energy dif- 
ferences between the paired leptons, i.e. (E 2 — Ei)(Ei — £3). In our notation, this product 
can be written, 

(E 2 - Ei)(E 4 - E 3 ) = 7172^1/92^1^2 cos 6»i cos 9 2 . (35) 

So this procedure places more importance on events with highly boosted Z bosons and/or 
events where the lepton is emitted along the line of the parent Z boson's direction of 
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Figure 15: The significances corresponding to the asymmetry Aq as a function of ^sm(b), 
for a Higgs boson of mass 150 GeV (left) and 200 GeV (right). We chose the other coupling 
coefficients a = 1 and c = 0. The inserts show the same quantities for a larger range of 
£sm(b). 

travel. Only the latter would affect our asymmetries and is similar (in principle and effect) 
to making a different choice of the observable, Oj. Thus, the procedure adopted by these 
authors is analogous to what we have done. 

In summary, using the observables 0i, O3, O4 and O5 and their corresponding asymme- 
tries at the LHC, we can in principle place limits on (or provide evidence for) CP-violation 
due to the simultaneous presence of CP-even and CP-odd form factors. For both a 150 GeV 
and 200 GeV Higgs boson, Oi, 3 , 4 and 5 all provide evidence (and 5 potential dis- 
covery) for some values of the additional non-SM couplings. Unfortunately, 2 and 6 are 
rather insensitive (due to the requirement of vector-axial interference) and cannot be used 
to place useful limits on additional couplings. The observables 3 and 5 can not unam- 
biguously rule out CP- violation, since their dependence also on 3?e(6*c) may provide an 
accidental cancellation with the terms proportional to a$t,e(c). However, 4 only depends 
on a9fte(c) and can thus test violation of the CP quantum numbers due to no n- vanishing 
a and c. With the three observables 03,4,5 at hand we can furthermore also extract the 
value of $le(b*c) and finally do consistency tests as well as reduce the effect of experimental 
errors. 

From a theoretical perspective, these asymmetries (if measurable) are sufficient to 
determine all the form factors a, b and c of our general CP-violating HZZ coupling, 
with real and imaginary parts. This is summarized in Table 1, which shows the various 
dependencies of the described observables on the form factors. We have 6 observables for 
the five unknowns a, 5fe(&), 5sm(b), 3fe(c) and £sm(c). If we furthermore assume that any 
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Asymmetry /form factor 


a 


Me(b) 


$>m(b) 


5fte(c) 


Qm(c) 


Ai 


X 
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A 2 
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(x) 


(x) 
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(x) 


A3 
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(x) 


(x) 
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(x) 


A A 


X 






X 




A 5 


X 


(x) 


(x) 


X 


(x) 


A 6 


X 




X 







Table 1: The dependence of the asymmetries Ai to Ae on the form factors a,b,c of the 
general HZZ coupling Eq. (2). (x) denotes a dependence which is suppressed if the 
additional form factors are small. 

product of b and c is very small (if b and c are loop suppressed or suppressed by some 
scale of new physics) we may neglect their simultaneous influence in the asymmetries A2, 
A3 and A$ and only require two of the asymmetries A3 to A$ (A2 not being of much use 
due to its smallness) to extract a and 9fe(c), while relying on A\ and A% for ^sm{c) and 
$sm(b) respectively. However, the analysis done here and the smallness of the asymmetries 
implies that only 9fe(c), and possibly £sm(c), are likely to be seen if non-zero. Also note 
that many of these asymmetries are highly correlated with one another. 

Of course, the final feasibility of detecting or placing limits on non-SM form factors de- 
pends on the real experimental environment. We simulated this here by taking the values 
given by the ATLAS studies. Any further refinement is beyond the scope of this study, but 
we have shown here the utility of these observables in providing unambiguous information 
on possible non-SM terms in the HZZ coupling and the consequent CP-violation. Any 
evidence or discovery of CP-violation crucially depends on the size of the non-SM form 
factors. Irrespective of this one may use these observables to place experimental limits on 
their values. 



5 Kinematical distributions as a probe of CP- 
violation 

As we have seen, the asymmetries discussed in section 4 are most useful for a Higgs 
boson with mass m# ~ 2m z- For a lighter Higgs boson the rates are much smaller and 
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the significance may not be sufficient for identifying CP-violation or setting satisfactory 
limits. In this case, one must rely on fitting shapes of kinematic distributions that depend 
on the CP character of the Higgs boson. From the discussions in the literature it is clear 
that the angle between the planes of the two fermion pairs coming from the decays of 
the Z bosons, and the polar angle of the fermions j\ or f 2 in the rest frames of the Z 
bosons, 6i (i = 1,2), are suitable variables [30] (see Fig. 1). 

1) The angular distribution in </> : In the decay process Eq. (1), let us consider the 
azimuthal angular distribution dT/d(j). Integrating Eq. (3) over 6i,6 2 and taking a CP- 
violating coupling with a and c non-zero 7 we find 

— — ~ bi + b 2 cos + 63 sin + 64 cos 20 + 65 sin 20 , (36) 
d<p 

where bi (i = 1, 5) are functions of m# and mz in terms of j a , 7&, 

h = a 2 (2 + 7 2 )+8|c|V 7 2 

= ~l$2 a mmia 
9tt 2 

63 = — aKe(c) 771772 X7a7fe (37) 
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2„,2„ ,2 



a 2 



64 = y-2| C |V 7b 2 

65 = — 2aKe(c) X'jb ■ 

Whereas the purely SM case (a = 1, b = c = 0) shows a distribution (see also Ref. [12]) 

dT 



d(j) 



1 + a 2 cos + a 4 cos 20 , 



° 2 = -^2 Vim 2T^ (38) 



1 1 



CZ4 = 



22 + 7 2 ' 

in the purely pseudoscalar case (a = b = 0, c 7^ 0) we have 

dr 1 

— ~l--cos20. 39 

d(f) 4 

In the CP violating case the inclusion of contributions from both the scalar and pseu- 
doscalar couplings alters the angular behaviour via the occurrence of sin0 and sin 20 
terms, and a reweighting of the other terms. Knowing the Higgs mass from previous mea- 
surements, any deviation from the predicted distribution in the purely scalar /pseudoscalar 
case will be indicative of CP violation. This can be inferred from Fig. 16 which shows the 



7 The expression with all three coupling coefficients a, b and c non-zero is given in the Appendix. 
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M H = 200GeV a= l,b = c = 0(SM) 

a = b = 0, c = i 

a=l, b = 0, c = i 




71/4 7C/2 371/4 71 

Figure 16: The normalized differential width for — > Z^Z — > (/i/iX/^) with respect 
to the azimuthal angle 0. The solid (black) curve shows the SM case (a = 1, 6 = c = 0) 
while the dashed (blue) curve is a pure CP-odd state (a = b = 0, c = i). The dot-dashed 
(red) curve and the dotted (green) curve are for states with CP violating couplings a = 1, 
6 = with c — i and c — i/2, respectively. 



azimuthal angular distribution for tjih = 200 GeV in the SM case, for a purely CP-odd 
Higgs boson and for two CP violating cases. The purely CP-odd curve will always show 
the same behaviour independently of the value of c since the curves are normalized to unit 
area. Therefore a special value of c could not fake the flattening of the curve appearing in 
the CP violating examples. This flattening even leads to an almost constant distribution 
in for the case c/a = i/2. It should be kept in mind though, that this method cannot be 
applied for very large Higgs masses where the dependence is washed out. One must also 
beware of degenerate Higgs bosons of opposite CP; since the decay products are the same, 
they will both contribute to the rate and must be summed coherently, possibly mimicking 
the effect seen above. 

This procedure is similar to that of Refs. [51,52] where log-likelihood functions were 
constructed and minimised to extract the coefficients in the vertex or yield exclusion 
contours. 

2) The angular distribution in 6i : Integrating Eq. (2) over and cos #2 provides a 
distribution in cos#i. For the CP violating case a,c^0, b = we find, 

dT 



a 2 [(t^ - l)sin 2 #i + 2] + 4|c| 2 x 2 7 2 (l + cos 2 #1) - 8a$m(c)r]ixj b cos 01.(40) 



d cos 9± 

In the purely SM case we recover, 

dT 



~ sin 2 ^-— ., (41) 

d COS 01 7a — 1 
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which for large Higgs boson masses (7 — > oo) reproduces the well-known behaviour 
sin 2 9\. In contrast, in the purely CP odd case we have 

dT 



1 + COS^ 01 . 



(42) 



d cos 9\ 

CP violation is manifest by a linear dependence on cos#i. However, due to the proportion- 
ality to rji the CP violating effect in the angular distribution is small, which is reflected 
also in the smallness of the asymmetry A\. See also the discussion in Section 4 and Fig. 3. 

2) The threshold distribution: In principle, information about the form factors of 
the HZZ vertex is also encoded in the dependence of its partial width on the virtuality 
of the Z-bosons [30]. In particular, looking at Eq. (9) one sees that only the term propor- 
tional to a 2 contains a linear dependence in (3. This is due to there being no momentum 
dependence in the SM HZZ vertex, in contrast to the additional non-SM terms of Eq. (2); 
the single f3 arises from the phase space. Consequently, one can distinguish a CP-even 
Higgs boson from a CP-odd Higgs boson decaying to ZZ* by examining the threshold 
behaviour since the CP-even excitation curve will be much steeper. This is illustrated in 
Fig. 17 where one can clearly see the steeper dependence on the virtuality M* of the most 
off-shell Z-boson for the CP-even case compared to the CP-odd case. 



0.04 



* 

Hi-i 



0.03 



0.02 



0.01 - 




30 35 40 45 50 55 60 65 
M* / GeV 



Figure 17: The normalized differential width for H — > Z^*>Z — > (/i/iX/^) with respect 
to the virtuality of the (most) off-shell Z-boson M*. The solid (black) curve shows the 
SM case (a = 1, b = c = 0) while the dashed (blue) curve is a pure CP-odd state (a = b = 0, 
c = i). The dot-dashed (red) curve is for states with CP violating couplings a = 1, b = 
with c = i. The vertical green line represents the nominal threshold at m# — mz- 



However, this behaviour near threshold will be dominated by whichever term has the 
lowest power of (3. So when one has a Higgs boson of mixed CP, the SM term will always 
dominate at threshold. This is also shown in Fig. 17 where the curve for the CP-violating 
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case sits almost on top of the SM curve near threshold. So while the threshold dependence 
is very good at distinguishing a pure CP-even Higgs boson from a pure CP-odd one, it is 
unfortunately not very helpful for distinguishing a CP-violating Higgs from the SM case. 

6 Conclusions 

In this work we have studied the process H — > ZZ^> — > 4Z, (I = e,/i) at the LHC to 
determine how well a general CP violating HZZ coupling can be tested. 

We examined the dependence of the partial width on non-SM form factors. By making 
use of the expected numbers of SM signal and background events, after cuts, provided 
by the ATLAS experiment, we produced exclusion plots for these non-SM form factors. 
We demonstrated that while large non-SM form factors may cause large deviations, it is 
difficult to distinguish their effect from an enhanced (or diminished) SM coupling. 

We then presented asymmetries which are non-vanishing when non-SM form factors 
are present in the HZZ coupling. We found a set of observables which, in principle, 
allows the extraction of the real and imaginary parts of all the complex form factors in 
the non-SM part of the HZZ vertex, if the significances are large enough. We analysed 
these asymmetries in the context of the ATLAS H — > ZZ^ — > 4/ study, and found that 
some of these asymmetries may be large enough to provide evidence of CP violation and in 
some cases even discovery, depending of course on the specific values of the CP violating 
contributions. In any case, these asymmetries will be useful in putting limits on any 
possible extra HZZ couplings beyond the tree-level SM, and deserve further experimental 
analysis. 

Furthermore, we presented an analytic formula for the partial width with full depen- 
dence on the final state azimuthal and polar angles, and demonstrated that the angular 
distributions may be exploited for Higgs boson masses below the threshold. Indeed, the 
azimuthal angle between the two decay planes of the Z bosons is sensitive to CP violation 
if the Higgs boson mass is not too large. 
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8 Appendix 

For the process H — > ZZ^ — > {fifi){f2h) with a general CP-violating coupling c.f. 
Eq. (2), we present here the differential distribution in the angle between the planes of 
the two fermion pairs coming from the decays of the Z^> bosons, taking into account the 
full dependence on the form factors a, b and c. The notation is as fixed in the text. 

(ir 

— ~ bi + b 2 cos + b 3 sin + 6 4 cos 20 + b 5 sin 20 , (43) 
a0 

where 

h = a 2 (2 + 1 2 a ) + \b\ 2 x 2 ^ + 8\c\ 2 x 2 1 2 + 2a^e(b)x lal 2 
9tt 2 

b2 = — ^T 7 ? 1 ^ 2 [a 2 l a + a^e(b) x%] 
9vr 2 

b 3 = -jg-771772 [^e(b*c) x 2 ^ + aKe(c) X'ja'jb] 
a 2 

U o I |2 2 2 

b± = — - 2 \c\ x 7 6 

6 5 = -2aKe(c) ^t, . (44) 

The polar angular distribution in 6\ is given by 
dT 



~ a 2 [( 7 2 - 1) sin 2 0! + 2] + |6| V 7b 4 sin 2 6 1 + 4|c| V 7 2 (l + cos 2 ^) 



d cos #1 

+ 2a3?e(6)a;7 a 7 2 sin 2 ^1 — 8a$5m(c)r]ia;7f, cos ^1 . (45) 
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